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We consider the first order differential equation with a sinusoidal nonlinearity and periodic time dependence, 
that is, the periodically driven overdamped pendulum. The problem is studied in the case that the explicit 
time-dependence has symmetries common to pure ac-driven systems. The only bifurcation that exists in the 
system is a degenerate pitchfork bifurcation, which describes an exchange of stability between two symmetric 
nonlinear modes. Using a type of Priifer transform to a pair of linear differential equations, we derive an 
approximate condition of the bifurcation. This approximation is in very good agreement with our numerical 
data. In particular, it works well in the limit of large drive amplitudes and low external frequencies. We 
demonstrate the usefulness of the theory applying it to the models of pure ac-driven semiconductor super- 
lattices and Josephson junctions. We show how the knowledge of bifurcations in the overdamped pendulum 
model can be utilized to describe effects of rectification and amplification of electric fields in these microstruc- 
tures. 



Pendulum and pendulum-like equations are ar- 
guably among the most important classes of equa- 
tions in modern nonlinear science 1 . The most 
often encountered representatives of this family 
may well be the driven and damped pendulum, 
6 + 7# + sm.6 = f(t), and its first order counter- 
part, the overdamped pendulum, jd + s'md = fit), 
the latter type being the topic of this paper. 
These equations appear, for instance, in the well- 
known Stewart-McCumber 2-3 and Aslamazov- 
Larkin 4 models of Josephson junctions. The si- 
nusoidal nonlinearity gives rise to a wide class of 
nonlinear phenomena that have important prac- 
tical applications: the ac-Josephson effect 5-6 and 
the modern voltage standard 7 are prime exam- 
ples of this. More recently, pendulum equa- 
tions have been found in the theory of semicon- 
ductor superlattices where they frequently oc- 
cur in the limiting cases of the governing dif- 
ferential equations^—. Moreover, overdamped 
pendulum equations are often encountered in 
mathematical models of synchronization of non- 
linear oscillators 15 . Further recent interest in 
the overdamped pendula has come from the field 
of high-T c superconductors: It has been demon- 
strated that stacked array of intrinsic Joseph- 
son junctions in magnetic field can be synchro- 
nized and described by overdamped pendulum- 
like dynamic o 16 i 17 . Properties of ac-driven over- 
damped pendula are also of importance in theo- 
ries of the amplification of microwave radiation in 
Josephson point contacts 1 ^—. And last but not 



least, in our previous work— we demonstrated 
how instabilities occurring in the overdamped 
pendulum are carried over to higher dimensional 
systems such as the strongly damped second or- 
der pendulum equation. Our present paper has 
two sides: mathematical and physical. Here we 
develop a mathematical technique which allows 
to find bifurcations in a class of overdamped pen- 
dula models for a wide range of their parame- 
ters, including a difficult but physically interest- 
ing case of low frequencies of driving force. We 
also show how this technique can be applied to 
symmetric physical systems demonstrating pen- 
dulum dynamics in some limiting cases. Our main 
focus is on the rectification and amplification of 
microwave radiation in unbiased semiconductor 
superlattices and Josephson junctions. 



I. INTRODUCTION 

We consider the first order ordinary differential equa- 
tion with a sinusoidal nonlinearity and arbitrary time 
dependence 

8(t) + G(t) sin 6 (t) =F(t). (1) 

The dynamics of the overdamped pendulum has been 
studied by several people, motivated by direct physical 
applications mentioned above. In spite of its apparent 
simplicity, novel nonlinear dynamics, e.g. strange non- 
chaotic attractors^ 3 - have been found. Here, we restrict 
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ourselves to a specific class of periodic forcing, specifi- 
cally consider bifurcations occuring in this systems, and 
its applications in physical systems. 

We will implicitly assume everywhere that F and G 
are real, continuous, and differentiable sufficiently many 
times. Our focus will be on functions F and G that have 
the following property 

F(t + T/2) = -F(t), G(t + T/2) = G(t). (2) 

With the above choice of external time-dependence, 
Eq. ([l} remains invariant under the transformation 

t^t + T/2, 9{t) -> -0(t + T/2) + 2fcvr, (3) 

and k is an integer. This type of forcing and the asso- 
ciated symmetry are of interest in many pure ac-driven 
physical systems, in particular bulk semiconductors and 
semiconductor superlattices, where breaking of symme- 
try ^ implies generation of a spontaneous dc bia a 8 ' 24 . 
In our previous work22 we in passing considered Eq. ([1]) 
with G(t) — 1 and F(t) — f coswi. We observed that the 
only instability that occurs is an exchange of stability be- 
tween two periodic solutions following symmetry and 
having the properties (9) = and (9) = w, where (•) 
stands for time-average across the period of the solution. 

We found that the symmetry breaking bifurcation, 
((9) ^ 0, 7r) 25 ' 26 of the strongly damped second order 
pendulum 

+ 70 + sin 6 = f cos art, (4) 

reduced to this instability in the limit of very large damp- 
ing. We conjectured that other ac-driven systems re- 
ducible to Eq. (P) undergo a type of bifurcation similar 
to the one found in the strongly damped pendulum equa- 
tion near the points where the exchange of stability oc- 
curs. This motivates our present extended study of the 
stability properties of system ([1} and applications to a 
number of physical systems. 

In this paper, our mathematical analysis is based on 
mapping of equation jl) to a particular second order 
linear differential equation. Such transformations have 
proved useful in the study of various linear and nonlin- 
ear differential equations. Following Priifer's application 
of the idea to Sturm-Liouville problems 2 ^, these changes 
of variables are sometimes called Priifer transforms. In a 
sense the reverse of this approach was taken in by Bon- 
deson et al£2- where the authors used a similar trans- 
formation to study quasiperiodically driven overdamped 
equation by relating the problem to a Schrodinger equa- 
tion with a quasiperiodic potential. We take essentially 
the same approach, but focus on the more specific prob- 
lem of periodically driven equation. 

On the other hand, applications considered in this pa- 
per are based on the connection of an exchange of sta- 
bility in the pendulum with the physical phenomena of 
amplification and rectification. Here we consider effects 
of microwave rectification and amplification in two pure 



ac-driven systems reducible to the overdamped pendu- 
lum: single-band lateral semiconductor superlattice and 
point-contact Josephson junction. 

Formally, by rectification we mean the conversion of 
pure ac excitation into response at even harmonics of 
some quantity that is a odd function of 9; for instance 
9 itself or sin 9. As an example of rectification, consider 
Eq. Q as a toy model where the drive f(t) corresponds 
to some ac applied field and that current is j(t) cx sm9. 
Rectification would then imply jdc — (j(t)) 0, i.e. 
obtaining a direct current response from a pure ac ex- 
citation, hence the term "rectification". For symmetric 
solutions (j(t)) = rectification is impossible, since 9 will 
only have odd harmonics, excluding possibly zeroth har- 
monic that is a multiple of n. Thus, symmetry breaking 
is a prerequisite for rectification. 

Note that rectification due to spontaneous breaking of 
symmetry in solutions, Eq. ([3]), should be distinguished 
from phase-dependent rectification due to breaking of 
symmetry in the equations 2 ^. The latter requires ex- 
plicitly introducing / that does not follow Eq. ©, for 
example an additional phase-shifted second harmonic 
cos(wi) -I- cos(2cj£ + 4>). 

For large damping we get overdamped first order pen- 
dulum for which exchange of stability arises for same 
parameters as symmetry breaking in second order pen- 
dulum. Rectification in pendulum is, however, rather 
artificial model which does not correspond to any real 
physical system. Nevertheless, in dynamical systems de- 
scribing realistic physical situations, symmetry breaking 
bifurcation is realized near values of parameters that are 
close to the values necessary for the exchange of stability 
in the overdamped pendulum. 

Here we consider a model of ac-driven lateral semicon- 
ductor superlattic o 11 ' 12 which is described by two first- 
order nonlinear balance equations which can be reduced 
to a sort of overdamped pendulum (see Eq. (|3H)) in the 
limit of strong nonlinearity. Symmetry breaking in bal- 
ance equations of lateral superlattice corresponds to rec- 
tification of applied ac electric fielcU^. We demonstrate 
how analysis of instabilities in Eq. (|TJ) can provide a quite 
useful information on the parameter space of rectification 
in these nanostructures. 

Our another application is related to amplification of 
infinitesimally weak signal in Josephson point contact de- 
scribed by Eq. ([1]), in which 9 - voltage, fit) - current 
and 9 itself is difference of phases of wave functions of su- 
perconductors in the junction. As a rule, an amplification 
of a small signal is observed near the onset of a dynami- 
cal instability- '- . Here we show a small-signal amplifi- 
cation near an exchange of stability and at frequencies of 
signal close to even harmonics of pump. Therefore, this 
effect of amplification can be considered also reminiscent 
of symmetry breaking bifurcation in strongly damped 
second order pendulum Eq. Q. Despite here even har- 
monics are forbidden by symmetry, nontrivial amplifica- 
tion of additional weak signal does exist at pump ampli- 
tudes and frequencies close to those necessary to realize 
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real symmetry breaking. 

The outline of this paper is as follows. In the next sec- 
tion we will introduce the change of variables that yields 
a second order linear differential equation and briefly re- 
capitulate on some known properties of its solution and 
their implications on Eq. (p}. We will then proceed to 
the more specific problem of forcing following Eq. ([2]) and 
show that an exchange of stability is the only instability 
occuring in this system. We then consider perturbations 
of the Eq. ([T]) and essentially prove our earlier conjec- 
ture that the exchange of stability is a limit of pitchfork 
bifurcations. In the subsequent section, we shift to a 
more practical approach: An approximate condition for 
the instability to occur will be derived in nontrivial case 
of large F and G. Finally, we go on to apply the re- 
sults to relevant physical problems. Technical details are 
presented in three Appendixes. 



II. EQUIVALENT LINEAR EQUATION 

We start by introducing a change of variables from 6(t) 
to new variables q\ (t) , <?2 (*) as 



6{t) = 2 arctan 



gift) 
«&(*) 



(5) 



We will denote the vector (q 1 (t),q2(t)) T by Q and the 
change of variables by 9{t) = C[Q(t)}. Since Eq. ((SJ alone 
does not fix the functions q\ , q% , we have some freedom in 
choosing the differential equations for the new variables. 
Here we opt for a particularly symmetric form of the 
equations 



d 
dl 



Q2 (i) 



1 / -G{t) F{t) 

2 I -F(t) G(t) 



<Zi(*) 
12 (t) 



(G) 



The coefficient matrix on the right-hand side of Eq. © 
will be denoted by A. We consider only periodic F and 
G, and therefore Floquet theory can be directly applied 
to the problem. We adopt the following notations for the 
Floquet solutions i = 1,2: 



= e B it 



Pi(t) 



= e Re{B i }* 



Pi(t), 



(7) 



where Bi are the (complex) characteristic exponents and 
Pi(t) are T-periodic functions. Functions Pi(t) contain 
the oscillating parts of the Floquet solutions. Additional 
lower indices will label the component of $j , Pi, and Pi , 
e.g. $i = ^i^) T ■ Due to the vanishing trace of 

A, the characteristic exponents have always the form (a) 
B\ 2 = ±-E>o + 27riki^/T, where Bq is real and are 
integers, or (b) B\^ — 2-Ki(±r + k\^)jT where r is real 
and not an integer. In the next section, we will show that 
the latter case is never realized if symmetry ([2]) applies, 
and thus case (b) will not be addressed in what follows. 
However, the special case of B\ t % — 2-Kiki^/T (Bq = 0) 
will be covered. 



The Floquet solutions give a complete description of 
dynamics of 9(t). Supposing case (a) from above holds, 
the general solution of the pendulum equation is 



9(t)=C 



cos — e 
2 



Pi 



. sin ^ e -Botp 2 



(8) 



where ipo, is a constant that depends on the initial value 
of 9, and whose value over — 7r < ip < n uniquely 
determines the solution 9 up to modulo 2tt. Clearly, 
when Bq ^ there are exactly two periodic solutions 9i, 
i = 1,2, that are simply given by the Floquet solutions: 



9i(t) = 2 arctan 



*i,2(t) 



(9) 



From Eq. © it follows that stable solutions of the over- 
damped pendulum correspond to the unstable solutions 
of the linear equation and vice-versa. This can be also 
seen from the relation that appears in Refs. HH and [32| 
and also applies here 



A 



1 

T 



G(t')cos9(t') dt' = 2\B 



(10) 



where A is the average exponential rate of growth for a 
infinitesimal perturbation of 9. Negative of its absolute 
value coincides with the maximal Lyapunov exponent of 
Eq. ([1]). Periodicity of the asymptotic solutions in the 
sense that 9(t + T) = 9 it) -\-2irn also follows immediately 
from the form of the Floquet solutions. 



III. SYMMETRIC CASE 

We now turn to our findings regarding Eq. (j6j) with 
forcing following Eq. ([2]). Matrix A transforms in the 
T/2-shift as A{t + T/2) = eA(t)e, where 



1 
-1 



(11) 



This property enables us to write the principal matrix U, 
U = AU, U (0) = /, on the latter half of a drive cycle in 
terms of the former 



U(t + T/2) = eU(t)eU(T/2). 



(12) 



The monodromy matrix M — U (T) can now be fac- 
tored into a square of the matrix M = eU(T/2), and the 
eigenvalue equation determining the Floquet solutions, 
M$i(0) = exp(SjT)$i(0), can be solved using the ma- 
trix M instead of M. Determinant of M equals —1, and 
thus the solution to the characteristic equation becomes 



j; arcsinh ( | tr M ] 
Bl.2 = { > ~ < 

I -fp arcsinh l| trM; 



(13) 



This shows B\ 2 are always real or real plus an integer 
multiple of 2iir/T. We will make frequent use of real part 
of B u Re-jBi} = B = 2arcsinh(tr M/2)/T. 
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For the periodic parts Pi of the Floquet solutions, the 
following now holds. Applying Eq. (fl"2"]l and M$j(0) = 
exp(B i T/2)$ i (0) one finds that P t {t + T/2) = eP;(t). 
Component-wise this property reads: 

P i>1 (t + T/2) = P iil (t), (14a) 
P i>2 (t + T/2) = -P ii2 (t). (14b) 

That is, the first component of Pi is T/2-periodic, and 
the second T/2-antiperiodic. It immediately follows that 
the periodic solutions 8i of Eq. ([1]) are symmetric in the 
sense of Eq. ([3]). Note that our definition for P>2, Eq. (p~3|) . 
includes an imaginary component, which contributes to 
the oscillating part of $2- Therefore, Pi is not real, and 
it is then more convenient to use the functions Pi instead. 
Now P\ — Pi, and so P\ has the same periodicity as P\ . 
On the other hand Pi = exp(2i7rf/T)P2, and thus P2,i is 
T/2-antiperiodic and p2,2 is T/2-periodic. 

Eq. (IT4|) also determines two properties regarding the 
rotations and the average value of the periodic solutions 
9i. Here we assume that roots of F are simple, that is, if 
F(t) = then F(t) ^ 0. We aim to connect the number 
of zeros of $i,2(i) over < t < T/2, here denoted n, to 
physically relevant properties of 9i - it will be shown that 
n indeed has significance to dynamics of real physical 
systems we are considering. Since 6i is symmetric, we 
can write 9i(T/2) — —9i(0) + 2irj, where j counts the 
positive direction crossings of the line 9 — n mod 27r. 
From Eq. © it can be seen that these crossings occur at 
simple zeros of <frj,2, and from Eq. (JXJ) that the direction 
of the crossing is given by the sign of F. Using again 
the symmetry, the average of 9i over t = . . .T, (9i), is 
jir. Now, it is easy to see that the parities of n and j are 
the same, and thus (9i) = nir mod 2tt. From Eq. (TT4")) 
it is clear that n is odd for 9\ and even for #2, and so 
(6>i) = 7r and (#2) = 0, both modulo 2-k. This shows 
that our previous finding^ regarding the averages of 9 
regarding the case F(t) — fs'mujt and G(t) = 1 holds in 
general. 

Further, n relates to the oscillations of sin^ and other 
quantities that are periodic 9. For instance, consider F is 
such that F(t) > (< 0) for < t < T/2 {T/2 <t<T). 
As 9 rotates it passes the upright vertical position n times 
and always in the positive direction, and so n gives the 
minimum and maximum number of oscillations of sin#j 
or cos#j in one half drive period. This has implications 
to the physical systems we are considering, since in these 
sin 8 and cos 9 have relevant physical interpretations. 

A. Exchange of stability 

Having established that the characteristic exponents 
are always real plus integer multiples of i27r/T, it then 
follows that the only possible type of instability is an ex- 
change of stability where one Floquet solution loses sta- 
bility and the other gains it. This in turn occurs when 
tr M — and as consequence B = 0. Noting that eigen- 
values of M are never equal, one finds that vectors $i (t) , 



i = 1,2, are linearly independent for all t and for any 
B . Thus, the Floquet solutions $i never map to same 
solution of Eq. (p} and the corresponding asymptotic so- 
lutions 9i(t) = C[$i(i)] never cross each other as a param- 
eter is varied. Consequently, the stability is exchanged 
without the solutions colliding, in contrast to a transcrit- 
ical bifurcation. We will later show, however, that the 
exchange of stability can be seen as a type of pitchfork 
bifurcation. From Eq. ([8]) it is clear that when Bq = 
all solutions to the overdamped pendulum are periodic. 
Since the superposition cos('0o/2)Pi + sin('0o/2)p2 that 
gives the general solution, Eq. , does not have period- 
icity analogous to Eq. (|14[) . the corresponding solutions 
9 are not symmetric. 

Further, crossing the instability has a clear effect on 
some relevant quantities. As was discussed above, the 
number of simple roots of ^.2 relates to the rotations 
and the average of 9. Let then $ + be the unstable Flo- 
quet solution, 9 + = C[$+] the stable periodic solution of 
the overdamped pendulum, and n the number of simple 
roots of over < t < T/2. Because at the instabil- 

ity switches between being <!>i and $2, n changes by 
one. Consequently, the average value of the stable peri- 
odic solution jumps by 7r, and further, since the value of 
n is intimately connected to the oscillations of sin 9 and 
cos 9, these quantities exhibit a change in in their fre- 
quency spectrum. We will later apply this finding in the 
section on lateral semiconductor superlattices. Clearly, 
the integer n partitions the parameter space into disjoint 
regions with the instability separating them. Thus, n 
serves as a convenient label for different regions of pa- 
rameter space. 

To help illustrate the bifurcation, we introduce a 
Poincare map n. Naturally we take this to be the stro- 
boscope map, defined so that 9(t + T) = H(9(t )) for 
some fixed <o which we take to be zero. If Bq ^ 0, the 
two fixed points 6*, i = 1, 2, of n are given in terms of the 
eigenvectors of M, or Floquet solutions at t = 0, <&;(()): 

"■*= 2 « c,a "(t!D- (i5) 

From the discussion above it follows that at Bq = 0, 
H(9) = 9 for all 9. In Fig. []] we have plotted a repre- 
sentative bifurcation diagram. We take G(i) = 1 and 
F(t) = /sinwi, where to = 0.3, and plot the fixed points 
9* as a functions of the forcing amplitude /. From the di- 
agram the bifurcation scenario can be easily visualized. 
At / = (outside the plot range) we have two fixed 
points: 9 = 0, tt, where the latter is naturally the un- 
stable point, since it corresponds to the upright position 
of the pendulum. At / ~ 1.356 we find the first bifur- 
cation. The initially unstable state becomes the stable 
one and vice-versa as the critical / is crossed. Exactly 
at the bifurcation, every point 9 is a fixed point of the 
Poincare map. From there on, increasing / further, we 
find the the fixed points exchange their stabilities again 
at / ~ 2.118 and / ~ 2.956, with the marginally sta- 
ble points spanning the whole phase space exactly at the 
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FIG. 1. Bifurcation diagram showing the fixed points of 
the stroboscopic map II of Eq. fTJ as a parameter is var- 
ied. Here G(t) = 1 and F(t) = /sinwt, where oj — 0.3. 
Forcing amplitude / is taken as the control parameter. Solid 
and dashed lines indicate stable and unstable fixed points re- 
spectively, while dotted line indicates marginally stable fixed 
points. These span the whole phase space and occur exactly 
at B = 0. 



bifurcation. 

The bifurcation described above bears resemblance to 
a pitchfork bifurcation. In fact it can be seen as a degen- 
erate pitchfork bifurcation (PB), since two new branches 
of fixed points emerge at the critical point. By degener- 
acy, we mean that these branches exist only exactly at 
the bifurcation point, span the whole phase space, and 
and are thus only marginally stable. This is in contrast 
to the (non-degenerate) PB, where the two additional 
branches of fixed points exist before or after the bifurca- 
tion, that are either stable (supercritical case) or unstable 
(subcritical case) . We note that this is reminiscent of the 
scenario observed for the second order, strongly damped 
pendulum^, where a (non-degenerate) pitchfork bifurca- 
tion was found near the criterion for exchange of stability 
in the overdamped equation. 

In terms of normal forms, non-degeneracy is under- 
stood as non-vanishing of a number of higher order 
derivatives with respect to the variable of the flow at 
an equilibrium^. We will see that this is indeed the case 
when we map Eq. (JTJ to an autonomous equation that 
is effectively a normal form on the circle . . . 27r. The 
form of the general solution, Eq. ((5J) suggests a natural 
way of mapping Eq. (JTJ into an equivalent autonomous 
form: the periodic parts P, describe the rotations of 6, 
while the exchange of stability is determined by the expo- 
nential factors and the superposition phase V>o- It seems 
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FIG. 2. Representative periodic solutions 8i, i — 1,2, of 
Eq. ((TJ with F(t) — fsinut, G(t) — 1 and the corresponding 
Floquet solutions $i just before [subfigures (a), (c), and (e)] 
and after a bifurcation [subfigures (b), (d), and (f)]. In both 
(a) and (b), the stable (unstable) periodic solution 9 is plotted 
with a solid (dashed) line. Similarly, Floquet solution map- 
ping to the stable (unstable) 9 is plotted with solid (dashed) 
line. Note that it is the unstable $ that corresponds to the 
stable 9. In the left-hand side subfigures (a, c, e), parameter 
/ = 2.9 which is just below the critical value of f — 2.956, 
while on the right-hand side (b, d, f) the parameter / is just 
above it, / = 3. 



natural to use a trial function where the constant phase 
"00 is replaced by a time-dependent tpit), that also ac- 
counts for the exponential factors: 



6(t) = 2 arctan 



rasfP M (i)+sinfP 2 , 2 (f) 



(16) 

We substitute this into Eq. ([]} and use Eq. ([6]) to obtain 
an equation for ip: 

ip = -2B Q miip. (17) 
No approximations were needed to derive Eq. (|17p . This 
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equation describes the approach to limit-cycles for the 
whole class of systems. Obviously, it is also the sim- 
plest non-trivial overdamped pendulum that has sym- 
metry @. Eq. (fTT|) has equilibria 0, n when B ^ 0, 
in the case Bq = every point is an equilibrium, cf. 
Fig. [TJ To compare bifurcations of Eq. (fl7|) to the pitch- 
fork bifurcation, we recall that normal form of a PB is 
r = r(/i ± r 2 ), where \i is the bifurcation parameter 33 . 
Prom Eq. (|17p . we see that near the equilibrium (7r), 
ip = ^ffit(,'(l — ip 2 /6 + •••)) where [i — 2B - Thus, we 
see that whereas in the non-degenerate case, only the 
leading order term vanishes at the bifurcation, here the 
right-hand side becomes identically zero at jj, = 0. 

To help visualize how the Floquet solutions relate to 
the solutions periodic solutions of the pendulum, we have 
plotted in Fig. [2] the periodic solutions 9i(t) — C[$j(f)], 
together with the Floquet solutions 3>i. We have used 
the same G and F as above and take / to be near the 
third bifurcation shown in Fig. [T] The left-hand side 
subfigures (a), (c), and (e) show the solutions for / = 2.9 
which is just smaller than the critical value of / = 2.956, 
while subfigures (b), (d), and (f) are plotted for / = 3.0. 
From (a) and (b) it can be seen that the Qi change lit- 
tle across the bifurcation point, only the stability is ex- 
changed. Similarly, the Floquet solutions remain roughly 
unchanged as the as the critical / is crossed, excluding 
the fact that the exponential envelope switches from de- 
caying to diverging or vice- versa. 



B. Effect of perturbations 

From Eq. (fT7|) it is evident that near a bifurcation 
the system is structurally unstable. An important ques- 
tion is then how dynamics change when a perturbation 
that allows for breaking of the symmetry ([3]), or explic- 
itly breaks is introduced. An exhaustive study of 
such perturbations is beyond the scope of this paper. 
Nonetheless we wish to show that the new type of dy- 
namics appear at the exchange of stability when a per- 
turbation is introduced. This is because to a large extent 
our motivation has been to show that the exchange of 
stability is in a sense a limit of bifurcations that occur in 
realistic physical systems that reduce to the overdamped 
pendulum. The perturbation is then to be understood 
as the terms removed from the original nonlinear sys- 
tem describing the real physical system to obtain the 
overdamped pendulum. Therefore we wish to show that 
symmetry breaking, or other type of dynamics appear 
exactly at the exchange of stability. 

We use trial function of the form given in Eq. (JTH]) 
to probe the response of the system to small additional 
terms. We introduce a perturbed system 

Q + G{t)wa.6 = F{t)+eH[9,t), (18) 

where < e < 1 and H{-9 + 2kir, t + T/2) = -H{6, t) 
does not necessarily hold. Substituting the trial function 



of Eq. (|16p into Eq. (TT8l) we obtain the equation 

= -2B sin V> - e£(V>, t) T ~{iP, t)H[6(iP, t), t], (19) 

where E(ip,t) = cos(ip/2)Pi(t) + sm(ip/2)P 2 (t). We have 
fixed the normalization of the Floquet solutions so that 
$i,i(0)$ 2 ,2(0) - $i, 2 (0)$2,i(0) = l 34 . Unlike in the case 
of Eq. (I16|) further approximations are needed. Smallness 
of e and Bq near a bifurcation can be used to simplify 
the above equation. 

First, we prove our earlier conjecture that the exchange 
of stability is the e — > limit of pitchfork bifurcations 
appearing in more realistic systems following Eq. (j2J or 
equivalent. We consider H = H(9) that follows the sym- 
metry H(—9+2kir) = —H(9), but contains even harmon- 
ics of 9. We note that the trial superposition S has the 
property E(-tp + 2kir, t + T/2) = (-l) k eE(ip, t). In this 
case, it then follows that also Eq. (fl9l) has symmetry ((2J. 
Using the averaging method 3 -^, we find that the averaged 
equation for ip will have the form 

ip ~ —2Bo sin ip + eai smip + ea-2. sm.2ip + ■ ■ ■ (20) 

where ak are constants which in general are nonzero. Co- 
sine harmonics of ip in the averaged equation are forbid- 
den by symmetry. Possible bifurcations of the perturbed 
system can then be qualitatively sketched by fixing ak 
and plotting the roots of ip = 0. As an example, if 02 7^ 
and ak = for k > 3, we find that the degeneracy of the 
pitchfork bifurcation has been lifted. A representative bi- 
furcation diagram is shown in Fig. El where the equilibria 
of Eq. ((201) are plotted for ea 2 = -0.1, a k = for k ^ 2. 
For comparison, the inset shows the degenerate limit of 
a-2 = 0. It can be seen that the degenerate PB is replaced 
by a pair of non-degenerate PBs that are supercritical for 
a 2 < and subcritical for a 2 > 0. Between the bifurca- 
tions, an equilibrium ip exists that is not equal to or n, 
and thus is a symmetry broken solution of Eq. (1201) . In 
this case, the t — > 00 solution 9 of the pendulum equa- 
tion will be described by a superposition of two Floquet 
solutions, and therefore, it too will not in general have 
symmetry ©. Thus we see that the exchange of stability 
is the limit of a symmetry breaking PB for symmetrically 
perturbed systems. 

If the perturbation H has only explicit, T2-antiperiodic 
time dependence, H = H(t) and H(t + T 2 /2) = -H(t), 
then the system will respond strongly when the perturba- 
tion frequency UJ2 — 2n /T a is close to an even multiple of 
lu = 2tt/T . This can be seen by noting that S T S has T/2- 
periodic zeroth and first cosine harmonic oiip- Therefore, 
the right-hand side of Eq. Q19[| will contain cosine har- 
monics of i\) whose coefficients oscillate at a frequency 
|2cj — CJ2I <C 1. Using again the averaging method, these 
terms will not vanish but contribute to slow, large am- 
plitude oscillations of ip. We will later show that this ef- 
fect has interesting consequences in the problem of weak 
signal amplification in Josephson junctions. Note also 
that this is in effect a dual of the pitchfork bifurcation 
described above - difference is that here the response fol- 
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-1.0 -0.5 0.0 0.5 1.0 

FIG. 3. Bifurcation diagram showing the equilibria of Eq. (|20[) 
for ea2 = —0.1, = for k 7^ 2. Horizontal axis is the bi- 
furcation parameter, [i = 2£>o- Stable (unstable) equilibria 
are plotted with solid (dashed) lines. Inset shows the degen- 
erate case of afe = for all k. Inset axis are the same as 
in the main figure, and the dotted line at fj, = shows the 
marginally stable equilibria. The main figure demonstrates 
that the degenerate PB is replaced by a pair of supercritical 
PBs at /j, — ±0.2. Between the bifurcations the equilibria 0, tt 
are unstable, while a new equilibrium with ip 7^ 0, n that cor- 
responds to symmetry broken solution of Eqs. ((20)) and ((T8|) . 



lows from a near even harmonic of the drive, not the 
angle 0. 

Finally, if the perturbation does not follow the symme- 
try @ but still depends on 0, one expects to see the co- 
sine terms appearing in Eq. (|20l) . Depending on the per- 
turbation there are several possible outcomes. The vari- 
ous bifurcation scenarios can then be enumerated by se- 
lecting the coefficients on the right-hand side of Eq. (l20l) . 
As an example, the pitchfork bifurcation may become 
an imperfect PB^, or essentially a saddle-node bifurca- 
tion, or the equilibria may be destroyed altogether near 
exchange of stability. 

In summary, we described the scenario for the devel- 
opment of only instability occuring in the system of type 
Eq. <[T|) following symmetry Eq. @. The scenario follows 
the one found in our previous work22. At a critical point 
a symmetric solution of Eq. ([T]) loses stability and an- 
other one gains it. The exchange of stability also marks 
the point where the average shifts from a minimum 
of the potential cos to a maximum or vice- versa. Ex- 
actly at the critical point neither of the Floquet solutions 
is diverging (and both indeed are still linearly indepen- 
dent), and hence the system will remain in some initial 
superposition indefinitely. In this special case, the sym- 



metry (j3|) need not be satisfied. This state is, however, 
only marginally stable, since a perturbation will neither 
decay nor diverge, and occurs only in a null set of pa- 
rameter values. 

We conclude by reviewing the some of the analytic 
methods of approaching the overdamped pendulum equa- 
tion to the exact linearization used here. The point of 
commonality to these methods is that they in one way 
or another consider the nonlinearity as a perturbation. 
For instance, the often used technique of using single har- 
monic trial function has requires the assumption that the 
sine term only has the effect of changing the amplitude 
and phase of the otherwise sinusoidal solution. Although 
the number of harmonics included in the truncation can 
be increased, the equations quickly get intractable. Av- 
eraging method can also be employed, as it indeed was 
in our analysis of Eq. (I20[) . however, again the poten- 
tial term needs to be small compared to the time scale 
at which varies. In contrast, the mapping to the lin- 
ear equation, Eq. (jB]) fully retains the nonlinearity whilst 
still allowing the use of tools for periodically forced linear 
equations, such as Floquet theory. 

Treating the nonlinearity nonperturbatively also allows 
us to also approach the limit of low frequency driving. 
In the following section, we derive approximate analyt- 
ical formulas for the solutions of the overdamped pen- 
dulum. Instead of studying the actual time dependence 
of 0, we continue with the focus on finding the critical 
points where the exchange of stability occurs. 

IV. ASYMPTOTIC SOLUTION OF EQ. @ 

We consider next approximate solutions of Eq. (TTJ) in 
the non-trivial limit of large F and G, or equivalently, 
low- frequency. We introduce a large parameter A into 
the problem by making the change F(t) — > \F(t), G(t) — > 
XG(t). In the leading order of A, we find that the equation 
we need to solve comes out as 

y+^(F(tf-G(t) 2 )y = 0, (21) 

where y = q\ or y = 52- Eq. (|21[) is found by taking 
the derivative of Eq. ([5]). Keeping only terms of order 
A 2 , one finds the equation (f2"Tj) for both qi and q2 sep- 
arately. We note that Eq. (|2"T|) docs not share the sym- 
metry of Eq. (JBJ), however, results of the previous section 
allow us to construct an approximate solution that has 
the expected properties. This follows from the fact that 
we need only to solve for the first half of a drive cycle 
t = ... Tj 2 and if need be, use Eq. (fl~2"j) to obtain the 
complete solution. 

Standard methods of asymptotic analysis^ can be ap- 
plied to Eq. (f2~Tj) to find its piecewise solution in the form 
(see Appendix |A1 . 

^(t)=|7^^(«^ 5l(i) +^-^ (t) )' ( 22 ) 
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where R(t) = {F{tf - G{tf)/A, &(*) = j£ \R{t)\ 1 ' 2 dt 

and w l = 1 (Wi = for R(t) < (R(t) > 0). We de- 

note by ti the turning-points (points such that R(U) = 0) 
and by N their number, i.e. i = 1, . . . ,N. Additionally, 
t = and tjv+i = T72. CoefBcients a,, 6, are solved from 
the initial conditions, and standard connection formulas 
for adjoining subintervals are applied. Using Eq. (f2"2"j) it is 
straightforward to construct a solution to any particular 
F, G. 

Although piecewise solutions naturally can be cum- 
bersome, we next show that tractable formulae can be 
obtained for quantities of interest. Naturally, we apply 
Eq. (|2"2"|) to calculating tr M, as its zeros define the crit- 
ical curves of the system. For simplicity, we limit the 
discussion to the case of two turning points, again, gen- 
eralizations are straightforward. With this restriction, 
the trace of M becomes 

tr M ~ -2sgnG(0) sinh((K + «i)A + In 2) cos ui A 
— 2 sgnG(O) cosh ((k — «i)A) sinwiA, (23) 

where n = £o(*i), i*>i = and m = &{T/2). This 

equation is one of the central results of this paper, as it 
allows for calculating the critical curves of overdampcd 
pendulum equations in the non-trivial limit of large G 
and F. 

In order to keep the treatment more concrete, we now 
fix G = 1 and F — fsmuit. We consider this restriction 
reasonable, since it was demonstrated in the previous sec- 
tion that differential equations of the form ([T]) exhibit 
the same structure as long as F, G have the appropriate 
symmetry. Thus, we can choose any such forcing as a 
representative of the class of equations we are consider- 
ing. Further, this choice has particular relevance to the 
physical systems we have in mind. 

Interestingly, with this choice, Eq. (|2ip becomes the 
Mathieu equation. In addition to constructions (|22|) 
and I23[) , we have the known solutions at our disposal. 
The trace of M comes out as 



(24) 



Here S = S(a, q; t) denotes the odd solution to the canon- 
ical form of the Mathieu equation^!, y + (a — 2qcos2t)y — 
0, with the (non-standard) initial condition S(0) = 1. 
The quantity tr M vanishes at parameters (a, q) for which 
S(t) is periodic, whether that period was n or 2tt. Val- 
ues of a corresponding to a periodic S(t) are the Mathieu 
characteristic values bk(q), k — 1,2,..., and thus the crit- 
ical curves are described by the equation 
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1,2,. 



(25) 



In addition to Eq. (|25| . which can be used to calculate 
the parameters for which tr M vanishes, for practical ap- 
plications we also need a way of computing tr M for any 
given (f,oj). Eq. is not well-suited for this purpose 




FIG. 4. Critical curves as functions of the parameters f,ui. 
Shaded regions indicate numerically found solutions of Eq. © 
[G(t) = l,F(t) — fsinuut] for which (6) — n. Critical curves 
given by Eq. (|25p are plotted with solid lines for ft = 1 ... 15. 
Dashed lines define the ten first Bessel roots [Jo(//oj) = 0], 
which are for clarity shown only for uj > / . 



since it requires the use of Mathieu functions with ar- 
bitrary parameters. Eq. (|23|) on the other hand has a 
more tractable form as it only requires the use of elliptic 
integrals: 

„ . , f^e{E{f 2 )} , \ Im{S(/ 2 )} 
tr M = -2 sinh — 1 u n + In 2 cos — 1 w n 



2 sin 



V u 
Im{£(/ 2 )} 



(26) 



where E is the complete elliptic integral of the second 
kind^. 

In Fig. @] we have plotted the critical curves as given by 
Eq. (|25l) together with the correct numerically obtained 
ones. For comparison, we have included a high-frequency 
approximation to the critical lines, Jo(//w) = where Jo 
is the Bessel J function of order zero22. The lines where 
f/u> equals a root of Jo are only in modest agreement 
for low frequencies, but improves as u) is increased. Our 
new result, Eq. (I25[) is, on the other hand, in very good 
agreement for low frequencies, and is accurate also in the 
opposite case of w ^ 1, especially for large /. Although 
not plotted, the condition trM = with M given by 
Eq. (|26[) also provides very good agreement to the com- 
puted critical curves. 

A quantity that will often be needed is the average of 
G cos 9, where 9 is the stable solution. Using Eqs. (|1Q[[T5|) 
we find the following equation that allows us to express 
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FIG. 5. Schematic figure showing the ac-driven lateral super- 
lattice. Electromagnetic wave is incident normal to the plane 
of 2D gas of conduction electrons. Electric field of the wave 
E(t) induces the current density j(t) flowing perpendicular 
to the superlattice layers. Along the direction of the current 
electrons experience a periodic potential with period o. 

(G cos 9) simply in terms of tr M as 

2 

(Gcos9) = 2IB0I = -|arcsinh(trM/2)|. (27) 

In summary, the main result of this section is Eq. (|23[) 
whose roots give the critical curves for F, G that are 
either large or depend slowly on time, obey symme- 
try and have exactly two distinct points ij^ta such 
that F(ti t 2) = ±G(t\ t 2)- For the case of constant G and 
sinusoidal F our three main results are (i) that inter- 
estingly the nonlinear equation reduces to the Mathieu 
equation, (ii) the Mathieu limit in turn allows us to write 
the critical curves of the system using the Mathieu char- 
acteristic values, with excellent agreement with the nu- 
merical data, (iii) Eq. (|2l>|) enables us to compute tr M, 
and consequently also (cos 9} [Eq. (|27|) ]. in the limit of 
slow external drive. This last result will be used in the 
following section. 

V. APPLICATIONS TO PHYSICAL SYSTEMS 

In this section we touch upon three physical prob- 
lems to which the theory developed above can be directly 
applied: Rectification of microwaves in lateral semicon- 
ductor superlattice o 11 ! 12 , amplification of high-frequency 
signals in Josephson point contacts, and modeling of 
Josephson junctions with critical current modulation^. 

A. Semiconductor superlattices 

We present the problem in the form it was introduced 
in Ref. H3 . A schematic figure showing the geometry 
of the problem is given in Fig. [SJ Plane electromag- 
netic wave incident on a lateral superlattice is considered. 
Electric field is polarized along the superlattice axis, that 
is, parallel to the direction of the current. The electron 
transport is studied by considering a single miniband 



with the standard tight-binding energy-quasimomentum 
dispersion relations^. The electron distribution follows 
the Boltzmann transport equation. From there, one is 
able to write ordinary differential equations for ensem- 
ble averaged electron velocity and energy. These form 
the well-known superlattice balance equation o 39 ! 40 . We 
present the equations here in their scaled form, in which 
the the maximum (minimum) electron velocity v and en- 
ergy w correspond to the value +1 (—1). Electric field 
inside the superlattice, denoted u, is also appropriately 
scaled. 

v = —uw — Tv, (28a) 
w = uv — T(w — w eq )- (28b) 

First equation of the set describes the balance between 
electron acceleration by the electric field and deceleration 
due to scattering, while the second describes the electron 
energy gain and dissipation due to scattering processes. 
The current density j is related to the average velocity v 
by j cx eN s v, where e is the elementary charge and N s 
is the areal density of 2D electron gas. The nonlinearity 
is controlled by the parameter T: T cx (7/iV^) 1 / 2 , where 
7 is rate of electron scattering - a high density of the 
electron gas corresponds to a large nonlinearity, or small 
r. Constant w cq is the scaled equilibrium energy, whose 
value we set to —1 for convenience. 

The interaction of the incident electromagnetic radia- 
tion with the conduction electrons is taken into account 
by employing the Maxwell equations with appropriate 
boundary conditions. Approximating the lateral super- 
lattice as an infinite conduction sheet, the scaled electric 
field entering the equations, u, becomes^ 

u = -u cos nt-T^v, (29) 

where uq and f2 are the amplitude and frequency of the 
external electric field. With the above, Eqs. (|28|) are 
rendered nonlinear. Unlike in the case of bulk superlat- 
tices where the relation between the total electric field 
u and the average velocity v is an additional differential 
equation^i, here the equation is algebraic. 

The overdamped pendulum equation is obtained via a 
formal change of variables, v = — A sin 9, w = — A cos 9. 
In the physically interesting limit of T <C 1 , the dimen- 
sionality of the system can be reduced (see Appendix IB"1). 
The following equation for 9 is obtained: 

g + f (cos9) F ) sm fl = u cosOt. (30) 
\ r (cos 9) ) 

One of the primary interests is the appearance of a 
spontaneous dc voltage. The term "rectification" was 
used for the conversion of applied ac irradiation into a 
dc field (u) ^ and current (j) cx (v) ^ 0, via the 
nonlinearity of these nanostructures. In terms of the 
variables A, 9, a prerequisite for such a current to ap- 
pear is that a limit-cycle does not follow Eq. Q, i.e. 
9(t + T/2) ^ -9{t) + 2kn. However, since the gov- 
erning equations are in fact symmetric in the sense of 
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Eq. ([3]), and based on the findings of Sec. HlIBl the break- 
ing of symmetry implies a pitchfork bifurcation. Conse- 
quently, rectification is expected in the real physical sys- 
tem when parameters are such that they correspond to 
the exchange of stability in the overdamped pendulum. 

In the earlier work 1 —, we have considered a simplified 
pendulum equations in which the contribution of (cos 9) 
was ignored, that is the overdamped pendulum Eq. ([1]) 
with G — 1 and F(t) = uq cos u>t. Our analytic analysis 
of that equation in|l2| was limited to the high driving fre- 
quency limit where the instability occurs in the vicinity of 
J (u /uj) = (cf. Fig. 2]). Comparing with results of nu- 
merical solutions of the superlattice balance equations, 
we observed that the rectification indeed exists nearby 
the Bessel roots. Here we apply the theory developed in 
the previous sections to Eq. (150"|) in order to find the re- 
gions of instability in a wider parameter space, including 
the case of low driving frequencies. 

We can solve the functional-differential equation (|30|) 
by considering the equation 9 + Ksin9 = uo cos fit. We 
require that K equals the coefficient of sine in Eq. ([50)1 : 



K 



(cos 6) (K) 



r 



(cos 8) (K) 



(31) 



Eq. (f2~T)) allows us to write (cos 9) in terms of tr M, and 
further, Eq. (l26l) gives trM using well-known special 
functions. Thus, roots of Eq. ([31]) can be easily com- 
puted numerically. 

Alternatively, parameter space structure of Eq. (|3"0")) 
can be studied by the following way. We introduce a 
simple change of variables (/, u) — > (uq,U): (uo,0) = 
{AT- 1 + A- X r) ■ (/», where A = (cos d) (/, w) = 
2\B (f,uj)\ can be found following Eq. (|2"7|) . Using this 
transformation, parameter space structure of Eq. (1301) in 
variables (uq, f2) can be studied for any T by calculating 
a single dataset of values (/, uj, B (f, uj)) from the pendu- 
lum equation ([1} with G = 1 and F(t) — f cos wi. For the 
case of T = 0.2 this procedure was applied to plot Fig. |6l 
where branches of the solutions with n = 1, . . . , 8 are dis- 
played. The pendulum limit of the superlattice balance 
equations [Eq. (|3U)) ] becomes invalid as (cos 9) < T, and 
thus we have chosen not to plot regions corresponding 
to (cos 9) < T. In the case of the simple overdamped 
pendulum, the position of, say, the nth parameter space 
region of solutions is determined solely by the external 
drive amplitude and frequency [Fig. [4]- In Eq. ([30| the 
average cosine affects values of (uo,f2) that admit a so- 
lution corresponding to a particular n. Account of the 
(cos 9) dependence has the effect of shifting the parame- 
ter space points corresponding to high values of (cos 9) to 
the higher values of u, which can be seen from the change 
of variables (/, uS) — > (uo,fi). For sufficiently low T, this 
shift is enough to make several regions with different n 
co-existing at fixed (uq,CI). The overlap of the different 
regions of solutions can be seen in Fig. [5] 

Comparing to Fig. 2 of Ref. 12|, we find that the 
branches n = 2,3,4 match the symmetry broken regions 




n 



FIG. 6. Parameter space of Eq. (|30p found for V — 0.2 employ- 
ing the asymptotic solutions discussed in Sec. IIVI Shading is 
used to distinguish different branches of solutions; white areas 
correspond to no valid solution ({cos 9} < V) or branch other 
than 1, .... 8. 



found using direct numerical simulation of the superlat- 
tice balance equations^. Branch n = 1 in Fig. [5] falls 
under the branch n = (not plotted), which suggests 
that a symmetry broken region might have been missed 
in the earlier work. Indeed, our subsequent numerical 
simulations confirm that there is a region of symmetry 
broken solutions associated with the instability between 
regions n = and n — 1 . It apparently was not detected 
because the initial conditions preferred the n = solu- 
tion. This immediately demonstrates the usefulness of 
our analytic results. 

The abstract theory developed in Sec. [ni] gives addi- 
tional insight into the dynamics of semiconductor super- 
lattice electrons. In Sec. IIIII we showed that the integer 
n was connected to the number of rotations of 9, and 
consequently oscillations of periodic functions depending 
on 9: in short, n counts the number of full oscillations of 
s'm8 and cos 9 across half the period of the drive. Since 
here sin 9 and cos 9 correspond to the average electron ve- 
locity and energy, respectively, we can now connect the 
dynamics of the charge carriers in the periodic poten- 
tial to the stability of the system. At the instability the 
integer n changes by one, and so the instability in fact 
marks the region where the state with n + 1 oscillations 
of energy or current become favourable to the state with 
n oscillations. Further, each of the branches n = 0, 1, . . . 
are characterized by the number of current oscillations 
across half of the drive period. 

In summary, the theory developed in the previous sec- 
tions devoted to the mathematical analysis of Eq. (Q] al- 
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lowed us to analytically probe the dynamics of Eqs. (|28|) 
in the physically interesting limit of T <C 1 , cj <C 1 • This 
limit is equivalent to very high nonlinearity, and has been 
largely inaccessibly analytically. Importantly, the analyt- 
ical results revealed the significant multistability in this 
system. We found that several branches of solutions can 
coexist at same parameters (uq,u)). The analytical re- 
sults also directly suggested, that a large region of recti- 
fication was missed. The region of rectification was found 
to be much larger than expected in the case of T = 0.2. 
Our new analytical findings are also consistent with the 
previous result that T < 0.4 is required for observing 
rectification. 



B. Josephson junctions 

As a second application, we consider the problem of a 
high-frequency gain (negative absorption) in microwave 
irradiated Josephson point contacts. The correspond- 
ing motion equation for the Josephson phase difference 9 
is^SI 

6 + sin = f sin u)t-\- e cos Qt. (32) 

In addition to the driving current / sin cut, a probe cur- 
rent I p = e cos fit has been added. The probe amplitude 
e is assumed small and the frequency f2 is incommensu- 
rate to the drive u. 

Having introduced a weak probe current I p , we wish 
to find the power absorbed by the junction, and its de- 
pendence on the frequency of I p . The absorbed power is 
described by A = (U ■ I p ), where U is the voltage across 
the junction. In terms of pendulum variables the absorp- 
tion takes the form Ajj = s(9(t) ■ sinftt). For Ajj < 
we have gain, i.e. a weak signal e sinilt will be amplified. 
The system placed in a cavity will radiate at frequencies 
for which Ajj < 0. Linearizing Eq. (13"2l and using the 
approach similar to PToL we find 

e 2 °° £1 2 
A " = J JL^ b ~* dk fi2 + ({cos e) + 2iku})2 > (33) 

where 

e ± / o *[cos0(t')-(cos0)] At' _ Q 2ikuit I Ok 

fc— — oo 

Derivation of Eq. (|3"3"|) is presented in Appendix [U] 

Near exchange of stability, (cos 0) — 2\Bq\ — > 0, the 
expression for absorption consists of terms of the form 
l/(0 2 — 4fc 2 w 2 ). These diverge when Q — > 2ku> (fc is an 
integer), and thus we expect to see strong gain for probe 
frequencies near an even multiple of the pump frequency 
when Bq ~ 0. 

Again, results of the previous section can be applied to 
find (/, uj) that result in such strong response. The range 
of parameters where gain occurs is rather narrow, and 
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FIG. 7. Regions of negative absorption computed from 
Eq. (|33[) . Shaded regions correspond to negative absorption 
for some probe frequency in the range fl = . . . 11.5 ■ w. 
Solid lines indicate exchange of stability as given by Eq. (|26|1 • 
Points labeled (a) and (b) indicate parameters for the absorp- 
tion profiles in the inset. Inset: Typical absorption profiles 
for low pump frequencies. Solid and dashed curves correspond 
to points (a), / = 1.5, u = 0.149, (b), / = 1.57, u = 0.13, 
respectively. 



hence an analytic first approximation is useful. In Fig. [7] 
we plot the regions of negative absorption using Eq. (|33[) 
together with the critical curves as given by Eq. (|26|) . 
Coefficients bk, dk were computed numerically from the 
Floquet solutions using the fact that the expansions in 
Eq. ([34]) are equal to (aP T P) ±1 [cf. App.[C], where Pis 
oscillating part of the unstable Floquet solution and a — 
(P(0) T P(0)) _1 . We also made several runs computing 
A3 j directly from Eq. (|3"2")l . While this direct procedure 
is in general more slow, its results are in a good agreement 
with the results obtained employing Eqs. (1331) and (IM)) . 

To illustrate the strong dependence of the absorption 
on the parameters we have also plotted typical absorp- 
tion profiles, dependence of Ajj on Q, in the inset of 
Fig. [71 the profile forms strong peaks just as the criti- 
cal curve for exchange of stability is approached. Gain 
peaks are located in the vicinities of even harmonics 2koj 
of the driving current. Therefore the system can radiate 
at frequencies close to 2 few. This effect is very different 
from ordinary generation of harmonics. Really, follow- 
ing symmetry ([3J overdamped pendulum can generate 
only odd harmonics (2fc + l)ui of the strong pump while 
even harmonics are forbidden. From physical point of 
view, generation at (exactly) odd harmonics is a spon- 
taneous process, whereas radiation at frequencies nearby 
even harmonics is stimulated emission. From the view- 



12 



point of bifurcations, the effect of gain in overdamped 
pendulum represents a reminiscent of symmetry breaking 
bifurcation in strongly damped second order pendulum 
Eq. (Q}, where the bifurcation results in appearance of 
even harmonics. Despite here even harmonics are forbid- 
den by symmetry, nontrivial amplification of additional 
weak signal does exist at pump amplitudes and frequen- 
cies close to those necessary to realize real symmetry 
breaking. 

Interestingly, similar gain profiles, centered near a 
characteristic Bloch or cyclotron frequency and its har- 
monics, attract much attention in physics of semicon- 
ductor superlattices, where them termed dispersive gain 
profile o 43 ! 44 . Moreover, the dispersive gain profiles were 
found in the models of ac-driven semiconductor superlat- 
tices as wel l 45 ' 46 . 

Recently, relatively strong coherent electromagnetic 
radiation of technologically important terahertz fre- 
quency band has been observed from dc-driven high- 
T c superconductors 4 ^ and arrays of niobium point 
contacts 4 ^. Physical mechanisms responsible for the ob- 
served radiation are still subject of intensive debates. 
Nevertheless, we notice that our simple model of a pure 
ac-driven junction gives a natural framework to general- 
ization of dc-drive configurations considered experimen- 
tally, and in this respect a further development of this 
theory may have a beneficial applied aspect. 

Finally, as a separate remark, we note that the general 
form of Eq. (fTJ) with G(t) ^ constant and the associated 
symmetry do have relevance to more complex models of 
Josephson junctions. Recently, a pendulum equation was 
used to model observed current-voltage characteristics of 
a type of Josephson junction with a strong applied mi- 
crowave field 3 ^. The magnetic field induced changes to 
the critical current were considered significant enough to 
be incorporated into the model. The resulting equation 
has, in the low-frequency drive case, a form correspond- 
ing to Eq. (Q]) and, in the absence of a dc current com- 
ponent, also the symmetry considered in current paper. 



VI. CONCLUSIONS 

We have studied the periodically driven overdamped 
equation ((T|) with periodic coefficients by using a type of 
Priifcr transformation. The linear form Eq. ^ allowed 
for easy analysis of the nonlinear system. We showed that 
if the system is driven by external forcing of the form 
Eq. ^ there exists only one type of instability in the 
system. The instability was identified as an exchange of 
stability between two periodic solutions, and it was found 
to be essentially a degenerate form of a pitchfork bifurca- 
tion. We showed that the degeneracy of the bifurcation is 
lifted when additional terms that also follow the symme- 
try are added, thus proving that the exchange of stability 
is in a way a precursor to symmetry breaking bifurca- 
tions in realistic physical systems. Also, we found that 
near the exchange of stability, the overdamped pendu- 



lum responds strongly to perturbations whose frequency 
is close to an even multiple of the drive frequency. We 
further used the linear form to find explicit solutions to 
the problem, and used them to construct a condition for 
the appearance of the instability. For the simple choice 
of F(t) = fsinut, G(t) = 1 the instabilities of Eq. §T§ 
can be well described by the Mathieu equation. We wish 
to emphasize, that up to the knowledge of the authors, 
no such low-frequency solutions have been constructed so 
far. 

From our purely mathematical findings there is then a 
natural cross-over into the field of physics. As was shown, 
near exchange of stability even weak perturbations can 
generate symmetry breaking bifurcations, or other types 
of strong response. Thus, regions close to this bifurcation 
in real physical systems are expected to show novel phe- 
nomena. We applied our methods directly the real phys- 
ical problems, namely rectification of pure ac irradiation 
in lateral semiconductor superlattices and amplification 
of a weak signal in Josephson point contacts. 

For lateral semiconductor superlattices, making use of 
our analytical methods, we were able to construct the 
parameter and phase space structure of this system. Ex- 
cellent agreement was found with our previous numeri- 
cal results 1 ^, with the additional finding that the system 
exhibits strong multistability: parameter space consists 
of overlapping regions each carrying a solution that is 
characterized by the extent of oscillations of the elec- 
tron gas average energy and velocity. Each of these re- 
gions was found to have an associated region of instability 
which in the low frequency drive and strong nonlinear- 
ity limit corresponds to rectification of the incident ac 
electric field - that is, the effect where the nonlinear in- 
teraction between the 2D electron gas and the external ac 
field generates a directed current and voltage across the 
superlattice. The analytical results revealed the complex 
way in which the multistability and the related regions 
of instability appear, and suggested that large regions of 
rectification were missed in earlier simulations. In fact, 
rectification persists for higher frequencies than was pre- 
viously thought. The new results supported the earlier 
finding that fairly high electron mobilities are required 
for rectification to appear. 

Second physical system that we considered was am- 
plification of a weak signal in Josephson point contacts. 
We found (i) a direct and striking correspondence be- 
tween gain, that is, negative absorption and the exchange 
of stability bifurcation - strong gain was shown to ap- 
pear in a narrow region just around the point where loss 
of stability occurs, when also the weak signal frequency 
was close to an even multiple of the strong ac driving 
current. Secondly, (ii) our analytical results make it pos- 
sible to find the regions of gain very accurately, or allow 
for computing them with little cost. Further, (iii) our 
findings on the dynamics of overdamped pendulum sug- 
gest insight into the physical process of amplification - 
the overdamped pendulum exhibits only odd harmonics 
which follows from the constraints set by the symme- 
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try, yet, we found that the system is expected to radiate 
when the weak signal to be amplified is tuned to even 
harmonics. 

Finally, we would like to point out intriguing similari- 
ties between gain profiles found here in the model of ac- 
driven Josephson junction (ac-driven overdamped pen- 
dulum) and dispersive gain profiles recently described in 
the models of ac-driven bulk [not lateral] semiconduc- 
tor superlattices i n 45 ' 46 . In these work o 45 ^ 46 , a superlat- 
tice is in essence modeled by two standard balance equa- 
tions for electron velocity and miniband energy 39 , that is, 
Eqs. (|28|) without v dependence in u. Taking into account 
some earlier findings^, we speculate that bifurcations 
of overdamped pendulum also determine conditions for 
high-frequency gain in ac-driven bulk superlattices, de- 
scribed by the standard balance equations, in the limits of 
high frequencies and rare collisions. A detailed compar- 
ison of properties of amplification in ac-driven Joseph- 
son junctions with the properties of dispersive gain in 
ac-driven semiconductor superlattice o 45 ' 46 in the limit of 
weak dissipation, however, goes beyond the scope of the 
present paper. 



proximately given by the usual WKB solutions 
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Appendix A: Construction of the asymptotic solution 



Here we briefly outline how the approximate asymp- 
totic solution of Eq. ((6|) is constructed using asymptotic 
solutions of Eq. (|21|) . Without loss of generality we can 
assume that R(t = 0) < 0. The interval t = 0...T/2 
is divided into subintervals Ij = (tj,tj+i) where tj, 
j = 1 . . . N is a turning point and N is the total num- 
ber of turning points. We additionally set = and 
tN+i — T/2. Solutions on each of the intervals are ap- 



V2k(t) 
V2k+l(t) 



(-R(t)) 1/4 

2fe+l . /. e m . f 



(Ala) 



A 



R(t) 
B 2 k 



R(t) 



^ cob (A&*+i (*) + £). (Alb) 



&(*)= / V\WY\dt'. 



(Ale) 



Each of the above solutions yk is only valid in its corre- 
sponding interval Ik - when R(t) < 0, solutions have the 
exponential form y 2 k and when R(t) > 0, the oscillatory 
form t/2fc+i is the appropriate solution [cf. harmonic os- 
cillator y + ry = 0, r real, solutions oscillate for r > 
and converge or diverge exponentially when r < 0] . 

Connection formulas for the coefficients Ck = 
(Ah, Bk) T can be derived by solving the problem at the 
turning points. For first order zeros of R(t), i.e. roots 
t* such that R(t*) ^ 0, the approximate solutions Xk 
around t = tk are given in terms of the Airy functions 
Ai, Bi as 

x k = A^]- 1 Ai(X 2 ^ k ) + B*SkY 1 ^{\ V Uk) (A2) 



with 



4>k 



-i?(t')] 1/2 dt' 



2/3 



(A3) 



where A^, Bj* are constants. Note that above we need 
to raise a complex number z to power 2/3, where z is 
either real or pure imaginary. Here, the argument of z 2 / 3 
is chosen to be or it, i.e. so that is real. Doing so 
it follows that if R(t) is increasing (decreasing) around 
tk, then 4>k{t), 4>k(tk) = 0, is continuous and decreasing 
(increasing) around tk ■ Away from the turning points the 
functions Xk asymptote into the WKB solutions given in 
Eq. (|A1[) . This allows one to write a linear relationship 
between Ck and Ck+i- C'k+i — WkCk, where 



W 2k = 
W 2k -i = 



2exp (K 2 feA) 




cosw 2 /c-iA -sincj2fc-iA 
sinw 2 fc-iA cosw 2 /c-iA 



1 ? w I ) (A4a) 
iexp(-ft 2 fcA) J 

(A4b) 



Here, w 2 fe-i = 6fe-i(*2fc) and K 2fc = ^kihk+i)- Denot- 
ing Y(t) — (y(t),y(t)) T , we wish to construct a matrix 
V so that Y(T/2) = V ■ Y(0). Above we have derived 
connection formulas for the superposition coefficients Ck 
across the whole interval 0...T/2, so all we need are 
matrices that map the initial values of Y , Y(0), to the 
coefficients Co, and the last coefficients Cn to the end val- 
ues Y(T/2). In other words, we need So and Si so that 
C = S ■ F(Q) and Y(T/2) = Si • C N . Using Eq. (jATall . 
(yo(0), yo(0)) T can be written in terms of the coefficients 
Cq = (Aq : B[)) t . After solving (Aq,Bq) the matrix So 
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can be read out. Similarly, using the same equation and 
writing out (yjy(T/2),yN(T/2)) T , the coefficient matrix 
Si is found. We obtain 



(A5a) 
(A5b) 




ri 








2Ar 








)■ 




















ro 




a 



Ar 
2 



n 
Mr- 



Appendix C: Derivation of formula for Josephson junction 
absorption 

We begin by linearizing Eq. ([52")) in e: Setting 9 - 
9 + 69, we get 



cos 959 — e cos fit, 



(CI) 



where 9 follows Eq. (|32j) with s = 0. The above has the 
exact solution 



where r = ^(O)! 1 / 4 , ri = R(0), and a = exp(Aftjv). 

Using Eqs. (TAT1) and ([A"3f the solution at t = T/2 can 
be obtained in terms of the initial values y(0),y(0). 



y(T/2) = ■ ■ ■ W S o Y(0). 



(A6) 



Finally, to get the value of Q at t = T/2, Eq. (fA"6|) is 
applied with the initial conditions 



y(0) = U(o)J' Y{0) - 

to get <7i(T/2) and qz(T/2), respectively. 



92(0) 
92(0) 



(A7) 



Appendix B: Overdamped pendulum limit of lateral 
superlattice balance equations 



In this appendix we present the derivation of Eq. p0[) . 
We start with the superlattice balance equations, 
Eqs. (|2"51) with an electric field following Eq. (|2"5)) . 



v = -(u in (t) -T- L v)w -Tv, (Bla) 
w = (iii n (t) — r _1 w)i> — T(w — w cq ). (Bib) 



Here, v, w are average electron velocity and energy scaled 
to dimensionless units and into range — 1 ... 1 , Ui n is the 
scaled incident electric field which is taken to be of the 
form Ui n = —uo cos fit. For details we refer the reader to 
Ref. [H. 

The pendulum form is obtained by making the substi- 
tutions v = — Asin9, w = — A cos 9. From Eqs. (|B1[) we 
get differential equations for A, 9: 



A = -FA - w cq T cos 9, 
9 = - - Wcq jJ smt 



(B2a) 
(B2b) 



We are interested in the case r <C 1. It follows that, since 
A a I, we may consider A a slow variable. Provided 
that r < fl we can say that to a first approximation A is 
constant over one cycle of the drive field T, T — 2-k jfl. 
Taking the average of Eq. (|B2al) we get A ss — w eq (cos 9) . 
Finally, setting w cq — — 1 for simplicity and substituting 
A -> (cos 9) into Eq. (|B2b|) we get Eq. ([50]), 



/*co S 0(t') df f e /*' cose(t") dt" CQsm , dt ,_ 



ee 



We note that since the periodic solutions 9 are symmet- 
ric, 9{t + T/2) = -0{t) + 2kn, then cos6> is T/2-periodic. 
Thus, we rewrite the exponentials appearing in 89 as 
a Fourier series with only even harmonics of fl. This 
gives Eq. (f34|). Using these expansions, and denoting 
(cos 9) = (3 for brevity, 59 becomes: 

56» = £^6 fe e i2a,fci e^* f d k ,eP t j 2u,k ' t cosflt dt 

k,k' 



£ b n - k d k e 



(f3 + 2ikuj) cos fit + £1 sin fit 
O 2 + (/? + 2ikuj) 2 ' 



(C3) 



Next, we calculate Ajj- Since cosflt harmonic is not 
present in 9, Ajj becomes e(cos fit ■ 59). By partial in- 
tegration, we get yljj = eO(sinrii • 59). The averaging 
then simply picks from 59 the coefficient of e 2mujt sin fit, 
n = 0, divided by two. Upon inspection of Eq. (|C3[) . 
Eq. ([33]) follows. 

Finally, we note that the expressions in Eq. are 
equal to (aP T P) ±1 , where a = (P(0) T P(0))- 1 . This fol- 
lows from the fact that 

e /*G(f)co.fl(f)df = a 'Q(t) T Q(t), (C4) 

where Q is any solution to Eq. ©, = C[Q], and 
a' = (Q(0) T Q(0)) _1 . The above can be verified by a 
straight-forward calculation of the logarithmic derivative 
of Q T Q. Specifically, for the unstable Floquet solution 
$ = cxp(\Bo\t)P we then have 



e J*G(f)coafl(t') dt' = a ^\Bo\tp^Tp^ 

Using Eq. (f2~T)) the claim follows. 
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